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We study the deterministic diffusion coefficient of the two-dimensional periodic Lorentz gas as a 

function of the density of scatterers. Resuhs obtained from computer simulations are compared to 

the analytical approximation of Machta and Zwanzig [Phys. Rev. Lett. 50, 1959 (1983)] showing that 

their argument is only correct in the limit of high densities. We discuss how the Machta-Zwanzig 

^^ ' argument, which is based on treating diffusion as a Markovian hopping process on a lattice, can 

^\ I be corrected systematically by including microscopic correlations. We furthermore show that, on a 

0^ ■ fine scale, the diffusion coefficient is a non-trivial function of the density. We finally argue that, on 

a coarse scale and for lower densities, the diffusion coefficient exhibits a Boltzmann-like behavior, 

whereas for very high densities it crosses over to a regime which can be understood qualitatively by 

the Machta-Zwanzig approximation. 



o 



O 

(N 



fl 
>. 



KEY WORDS: deterministic diffusion; periodic Lorentz gas; computer simulations; random walk; 
chaotic scattering; symbolic dynamics; Boltzmann approximation; density expansion of transport 
coefficients 



> ■ PACS numbers; 05.40.-a,05.45.-a,05.60.-k,51.10.+y 

00 ■ 

m 

o 

'^ ; I. INTRODUCTION 

a^ 

^^ . One of the central themes in the theory of chaotic transport is the problem of deterministic diffusion. Over the 

past several years, deterministic diffusion coefficients have been computed for a variety of low-dimensional model 

systems. The most elementary ones are chains of one-dimensional maps [0-H, which have been studied with cycle 

'^ ' expansion methods P-@l . By applying techniques based on Markov partitions it has been found that the deterministic 

' ' diffusion coefficient of a simple pieeewise linear one-dimensional map is a fractal function of the control parameter 

rrt oty • Long-range dynamical correlations in the microscopic chaotic scattering processes of the system are responsible 

(— ! for this surprising behavior |10|-|l^ . More coniplicated models exhibiting deterministic diffusion are two-dimensional 



Q maps such as multi-Bakers |l^lq1, cat maps 17|, standard maps |_8|, and related models [19|]. Also in these models 

L^ the diffusion coefficients are fractal jlfl] , or are highly nontrivial functions of a control parameter p7|-p^ . 

• i-H . These simple maps share fundamental physical and mathematical properties of more complicated dynamical systems 

^^ ' such as the periodic Lorentz gas EG]. For this model, the diffusion coefficient has been computed with different 

H methods: in Refs. ||2l]-[2i^ it has been obtained via the Green-Kubo formula on the basis of computer simulations, in 

■ - - ' Refs. p^,pH| techniques of cycle and periodic orbit expansions have been applied, and in Refs. fsl,^ the diffusion 

coefficient has been obtained from escape rate methods as well as from the Hausdorff dimension of the fractal repeller 

of the Lorentz gas. 

In this paper we focus on the behavior of the diffusion coefficient in the two-dimensional periodic Lorentz gas 
under variation of the density of scatterers. In Section II we introduce the model, briefly sketch the simple analytical 
approximation obtained by Machta and Zwanzig pl| and compare it to our numerical results obtained from computer 
simulations. In Section III we explain how the Machta-Zwanzig argument can be corrected systematically by taking 
microscopic correlations into account. In Section IV we compare our numerical diffusion coefficient to a simple 
Boltzmann approximation and argue that, on a coarse scale, there is a dynamical crossover for the diffusion coefficient 
as a function of the density. More detailed computer simulations show that the diffusion coefficient exhibits a non- 
trivial flne structure as a function of the density, as is discussed in Section V. In Section VI we summarize our results 
and relate them to analogous findings in one- and two-dimensional maps. 



II. DIFFUSION AS A SIMPLE MARKOV PROCESS: MACHTA-ZWANZIG APPROXIMATION AND 

NUMERICAL RESULTS 



The periodic Lorentz gas is a standard model in the field of chaos and transport l2flJ2q-pO|. It mimics classical 
diffusion in a crystal and is isomorphic to a periodic system of two hard disks per unit cell ||31[ . The geometry of 
the model is shown in Fig. y A point particle of mass m moves with constant velocity v in an array of circular hard 
scatterers of radius R arranged on a triangular lattice. Upon collisions with the scatterers the particle is reflected 
elastically. In the following we use units for which v = 1, m = 1, and R = 1. The lattice spacing of the disks is then 
2 + w, where w is the smallest inter disk distance. The gap size w is geometrically related to the number density n 
of the disks by 

n = 2/[V3{2 + wf] . (1) 

The gap size w, or the number density of the disks n, respectively, is the only control parameter and completely 
determines the dynamical properties of the system. At close packing, i.e., w = 0, the moving particle is trapped in a 
single triangular region formed between three disks (see Fig. |l|). For < w < Woo — 4-\/3 — 2 = 0.3094, the particle 
can move across the entire lattice, but it cannot move collision-free for an infinite time. For w > Woo the particle can 
move arbitrarily far between two collisions. Here, Woo denotes the gap size at which the particle first sees such an 
"infinite horizon" . Bunimovich and Sinai proved that for < w < Woo the system is ergodic, and that the diffusion 
coefficient exists, while it diverges for w > Woc [ p2||3^ . 

In Ref. |£l|] Machta and Zwanzig have derived a simple analytical approximation for the diffusion coefficient D. The 
basic idea is that diffusion can be treated as a Markovian hopping process between the triangular trapping regions 
indicated in Fig. |l|. For this purpose, they have calculated the average rate r~^ at which a particle leaves such a trap. 
According to a simple phase space argument, this rate is determined by the fraction of phase space available for leaving 
the trap divided by the total phase space volume of the trap. Furthermore, for random walks on two-dimensional 
isotropic lattices the diffusion coefficient is D = P/{At), where I = {2 + w)/\/3 is the distance between the centers of 
the traps. This leads to the Machta-Zwanzig random walk approximation of the diffusion coefficient 

^^^-n[V3i2 + wy~2n] ' ^'^ 

This equation is shown as the dotted line in Fig. g. Included in this figure are single data points for the diffusion 
coefficient obtained by various authors using different methods: The filled circles are numerical results of Ref. pl| , 
the stars are from Ref. |^. These data points have been obtained by employing the Green-Kubo formula, where D is 
determined from an integral over the velocity autocorrelation function. The squares have been computed in Ref. pi] 
by periodic orbit expansions, the empty circles and triangles are from Ref. [p7|, where they have been computed by 
applying the escape rate formalism, and via the fractal dimension of the repeller of the Lorentz gas, respectively. The 
crosses connected with lines represent our new data which we have obtained from computer simulations by means of 
the Einstein formula 

Z,.li.<^ , (3, 

where the brackets indicate an ensemble average after time t. We calculated (Ar^(t)) by averaging over long trajecto- 
ries. Depending on the density n each trajectory contained from 6 x 10^ to 3 x 10^ collisions. After a short transient 
(Ar^(i)) grows linearly with a slope of AD. We determined the diffusion coefficient D by fitting a straight line to 
(Ar^(t)) in the linear regime. Though in principle equivalent with the Green-Kubo formalism, the Einstein approach 
is numerically more efficient in the Lorentz gas, where the equations of motion can be integrated exactly. 

Fig. shows that the Machta-Zwanzig approximation Eq. (0) is only correct in the limit of small gap sizes w — > 0. 
For larger values of w and up to approximately w < 0.1 Eq. (0) clearly overestimates the exact diffusion coefficient, 
whereas for w > 0.1 it significantly underestimates diffusion. In the following section, we will discuss the physical 
reason for this failure of the Machta-Zwanzig approximation and how it can be corrected. Note also that for w -^ Woo 
there is no evidence for any singularity in the diffusion coefficient reminiscent of critical behavior. In Ref. |g3[, the 
same observation has been made based on the analysis of velocity autocorrelation functions. For explaining this lack 
of criticality the authors have argued that this transition rather represents a simple geometric feature of the system 
than being a consequence of long-range correlations. 



III. CORRECTION OF THE MACHTA-ZWANZIG APPROXIMATION: CORRELATED MICROSCOPIC 

SCATTERING 

A. Collisionless flights across a trap 

As Machta and Zwanzig remarked themselves, for larger w there is a non- vanishing probability pd for the particle 
to move collision- free across a trap. In Fig. y (a), we have calculated this probability from computer simulations as 
well as in a straightforward analytic al a pproximation, which is described in Appendix Kl Note that collisionless flights 



occur only for w > 0.1547, see Eq. ( |A6D . 

If we rely on the Machta-Zwanzig picture of diffusion as a hopping process with frequency t~^ over distances I, we 
can correct Eq. (0) in the following way: If aparticle moves collision-free across a trap, it travels within the time r 
over a larger distance than assumed in Eq. (||). For this larger distance we take the distance between a center of a 
trap and the center of its next nearest neighbor which is I2 = v3^- These processes would thus yield a larger diffusion 
coefBcient Di^ = Z|/(4r) = SDmz, where Dmz is the diffusion coefficient of Eq. (||). We now define the corrected 
Machta-Zwanzig diffusion coefficient Dd by weighting the contribution of collisionless flights via the probability p^ 
of Fig. I (a), 

Del == [1 -Pcf]DMZ +Pc[iDMZ 

= [1 + 2pci]DMZ ■ (4) 

The result is plotted in Fig. pi (b). For w > 0.1547, the revised formula Eq. (0) improves the original Machta-Zwanzig 
approximation considerably, however, it is still much smaller than the numerically exact results. 

B. Probability of backscattering 

A second significant contribution to the correction of Eq. (||) is determined by the backscattering probability pbsj 
which is the probability of the moving particle to leave the trap through the same gap where it entered. We computed 
Phs numerically by repeatedly injecting the particle through a specific gap and observing trough which gap it left the 
trap. The particles are initially situated at one entrance of a trap, and they are uniformly distributed in the respective 
phase space of this entrance, which consists of the position of a particle on the entrance line, — 1 > 2x/w > 1, and 
of the sine of the angle between the velocity direction and an axis perpendicular to the entrance, — 1 > sin a > 1. 
Alternatively, pbs can be determined from a single long trajectory. The backscattering probability obtained from our 
simulations is shown in Fig. H (a). 

The Markovian approximation of Machta and Zwanzig Eq. (Q) corresponds to a backscattering probability of 1/3. 
However, for small values of w the numerically computed probability is significantly larger than 1/3, whereas for larger 
w it is considerably smaller. The reduced probability for backward scattering at larger w is in part due to collisionless 
fiights across the trap. It is remarkable that the backscattering probability is different from 1/3 even for gap sizes 
where the average number of collisions between inter-trap hops is large. At w ~ 0.055, where the backscattering 
probability reaches its maximum value of phs — 0.38, more than 17 collisions occur before the particle hops to the 
next trap and for w = 0.02, where pbs is still close to 0.36, the number of collisions is greater than 50. Correspondingly, 
the detailed functional form of pbs appears to be quite intricate: Below the maximum at w = 0.055 there are at least 
three regions where phs decreases approximately linear in w with different values of the slope. As the numerical results 
indicate, the function must eventually drop extremely sharply to pbs(O) = 1/3. Details of these regions are shown in 
the magnification included in the figure. The vertical lines separating different regions correspond to the respective 
lines separating regions of different slope in the main figure. A very close look reveals that the fine structure of all 
these different regions appears to be quite similar, however, we note that this structure is essentially within the range 
of our numerical errors. 

Fig. |5| shows the initial conditions in the (a;,sina)-plane leading to backscattering for a gap size of ■«; = 0.1. Here, 
the X-axis is parallel to the trap entrance and the origin is located at the center of the entrance. Initial conditions 
were drawn at random from a uniform distribution in this plane and a dot was plotted whenever the particle left 
the trap through the same aperture where it entered. Simulations at different densities yield similar pictures. The 
complicated, intertwined structure is reminiscent of a fractal set. The modifications of this structure by varying w 
must be related to the changes seen in the function pbs of Fig. (a). This may explain why, even on the coarse scale 
of H (a), Phs is not a simple function of w. Moreover, the detailed changes of this fractal set may be reflected in the 
respective detailed changes of pbs on the fine scale, as illustrated in the magnification. 



Having the probability of backscattering pbs we can now perform a second correction of the original Machta-Zwanzig 
diffusion coefficient Eq. ^). For this purpose, we again assume that diffusion can be treated as a hopping process 
with a frequency t~^ over distances I. We now inquire to which traps the particle can move by performing two jumps 
within a total time interval of 2t. There are only two possibilities: either the particle suffers backscattering, that 
is, it goes back to its original trap and does not contribute to any actual displacement within 2t, or it moves over a 
distance I2 to the left or to the right of its original trap. Thus, the corresponding diffusion coefficient reads 

[1-Pbs]^i 

= [1 - Phs]3/2DMz . (5) 

-Dbs is plotted in comparison to our Einstein formula results, and to the original Machta-Zwanzig Eq. (0), in Fig. 
(b). For smaller w, the corrected diffusion coefficient approximates the numerically exact values quite well. Thus, we 
conclude that the existence of backscattering is basically responsible for the Machta-Zwanzig argument overestimating 
diffusion for small values of w. For larger w Eq. (^ again improves the original Machta-Zwanzig approximation, like 
the previous approximation Eq. (0), however, like this it yields much smaller results than the correct values. 

C. Combining collisionless flights, backscattering, and a symbolic dynamics 

So far we have identified two different microscopic scattering mechanisms which were not included in the original 
Machta-Zwanzig approximation for the diffusion coefficient, and which evidently play a significant role for under- 
standing the full density-dependent diffusion coefficient: Namely, the probability of collisionless flights across a trap 
leading to the diffusion coefficient of Eq. (Q), and the probability of backscattering leading to Eq. (^. Consequently, 
one may think of combining these two different approximations within a single diffusion coefficient formula. This can 
be performed by simply replacing the Machta-Zwanzig diffusion coefficient Dmz hi Eq. (|^) by Dd of Eq. (Q) yielding 

Di = 3/2[f - pbs] [1 + 2pci]Dmz ■ (6) 

This function, denoted as a first order approximation, is shown in Fig. M (a) in comparison to the numerically exact 
results and to the Machta-Zwanzig Eq. (g). For larger w, this combined approximation is much closer to the numerical 
results than the original Machta-Zwanzig formulation, however, there still remains a notable quantitative difference. 

In the same figure, corrections of higher order based on the idea of Eq. (|g) are depicted. In the following we just 
outline the basic concept, the detailed formulas are then given in Appendix M. The higher-order corrections of the 
diffusion coefficient are obtained by numerically computing the probabilities of higher-order backscattering, and by 
building them into a respectively generalized version of Eq. (ph. This generalized expression of the backscattering 
diffusion coefficient I?bs is then combined with the respective collisionless flight-diffusion coefficient Dd of Eq. (^ in 
the same way as before. The probabilities of higher-order backscattering have been computed on the basis of a simple 
symbolic dynamics, as it can be defined in case of simple backscattering: We followed a long trajectory of a particle 
in the Lorentz gas. For each visited trap we labeled the entrances through which the particle entered with z, the 
exit to the left of this entrance with I, and the one to the right with r. Thus, a trajectory in the Lorentz gas can be 
mapped to a sequence of symbols z, I, and r. Note that the symbolic dynamics we are using here is different to the 
one applied in other work in that we are labeling the three gaps of a trap, whereas in previous work the single disks 
accessible after a collision have been chosen, which required an alphabet of 12 symbols [p|j2C|]. p{z) = pbs is then the 
backscattering probability depicted in Fig. (a), whereas, due to symmetry, p{l) — p{r) = (1 — p{z))/2 corresponds 
to forward scattering. Combining these probabilities with the correction by collisionless fiights yields Eq. (||) as a first 
order approximation. 

Higher-order correlations can be calculated by taking into account the probabilities of longer symbol sequences. 
For example, the second order approximation involves probabilities corresponding to nine symbol sequences each 
consisting of two symbols, p{zz), p{zl), p{zr), p{lz), p{ll), p{lr), p{rz), p{rl), p(rr). In complete analogy to Eq. (|^), 
a respective second-order approximation of the diffusion coefficient is then computed on the basis of these numerical 



probabilities, and by associating to them the respective distances traveled within a time interval of 3t, see Eq. (Bl). 
In the third-order approximation the probabilities correspond to sequences of three symbols, for example, three times 
back scattering within a time interval of 4r corresponding to p{zzz), etc., and lead to the diffusion coefficient of Eq. 
(B2). In general, the fc-th order approximation involves the probability of sequences of k symbols. Fig.|| (a) shows 
that by including such correlations the respective higher-order approximations of the diffusion coefficient converge, 
and basically move closer to the numerical exact results. However, for u; > 0.1 the approximations converge to 
the exact results apparently only globally, and not locally. A clear sign of this is that in third order the respective 



approximation already exceeds the numerical data for large w, although the approximations seem to approach the 
exact data generically from below. On the other hand, one should take into account that this approximation scheme 
is based on a purely heuristical ansatz for the diffusion coefficient. Thus, an exact convergence should not necessarily 
be expected. 

We therefore studied the effect of increased (or decreased) backscattering with the help of a lattice gas computer 
simulation. In such a simulation the Lorentz gas is mapped to a honeycomb lattice where the sites of the lattice 
represent the traps. The moving particle hops from site to site with frequency t~^ , which is identical to the exact 
hopping frequency used in Machta-Zwanzig theory. We first describe the first order approximation: At each step the 
particle hops back to the site where it came from with probability phs or to one of the other sites with probability 
(1 — pbs)/2. The backscattering probability pbs used in the lattice gas simulations is the one numerically obtained 
from simulations in the Lorentz gas. Also in the lattice gas simulations we determine the diffusion coefficient from 
the Einstein formula Eq. (||). The results of these simulations are shown in Fig. H (b) by a dotted line. Taking into 
account these first order correlations brings the diffusion coefficient closer to the correct value, but there still is a 
considerable deviation. This indicates the importance of higher order correlations. 

These higher order correlations can be obtained by correlating more than two hops between neighboring traps. To 
determine the diffusion coefficient for such multiple hopping events by lattice gas simulations we used the probabilities 
p{lrz) calculated up to fourth order from long trajectories in the Lorentz gas. The results of these calculations are 
shown in Fig. @ (b). As can be seen in the figure, the diffusion coefficient obtained from this scheme converges very 
quickly to the numerically exact results, in particular for small w. For larger w, the convergence is somewhat slower, 
however, the fourth order approximation can be hardly distinguished from the numerically exact results on the scale 
of Fig. g (b). For growing order the diffusion coefficient converges to the correct results exactly. 

IV. BOLTZMANN APPROXIMATION FOR THE DIFFUSION COEFFICIENT 

In the previous section, we have discussed how the Machta-Zwanzig theory can be improved systematically. An 
alternative approach for understanding diffusion in the Lorentz gas may be obtained from kinetic theory, which 
essentially applies at low densities. In Refs. Q,^, the Boltzmann approximation for the random Lorentz gas, where 
the scatterers are distributed randomly in the plane without overlap, has been used to compute the diffusion coefficient 
to 

3 

DBo{n) = -lc{n)v , (7) 

o 

where Ic — l/{2Rn) is the collision length of the moving particle. At higher number densities n the excluded volume 
of the scatterers becomes important, and the respective shorter collision length reads 

DBoin) ^ ^{- - n) , (8) 

lb n 

where we have changed to units with v = 1 and R = I. The same result for the collision length can be obtained 
directly from the Machta-Zwanzig argument of Section II |2l). In case of the periodic Lorentz gas n is related to 
the gap size w by Eq. (m). Note that for w = 0, where the particle is trapped, the corresponding density is n < cx3 
and the collision length is still finite. Thus, Eq. (^) implies a diffusion coefficient of Dbo{w = 0) > 0, that is, this 
approximation inevitably leads to an offset between Dbo{w = 0) and the exact result of D{w = 0) = 0. Therefore, it 
can immediately be concluded that I?bo is not correct at the highest densities. 

In Fig. M, Eq. m) is shown in comparison to the numerically exact results and to the Machta-Zwanzig approximation 
Eq. (||). For this purpose, the diffusion coefficient has been plotted as a function of 1/n, and an offset of AD = 0.025 
has been subtracted from Eq. (g). According to this equation, -Dbo(?i) should go linear in 1/n with a slope of 3/16, 
and the offset AD is the only fit parameter. Indeed, Fig. ^ shows that, on a sufficiently coarse scale and over a wide 
range of lower densities, D matches quite well to this functional form. In particular, the value of the slope appears 
to be almost exact. Only for larger densities D deviates from linearity in 1/n. However, here the Machta-Zwanzig 
approximation seems to describe the functional form at least qualitatively quite well. 

From a conceptual point of view one may inquire why a Boltzmann description should at all be applicable to the 
periodic Lorentz gas. The basic difficulty is that because of the occurrence of an infinite horizon, a diffusion coefficient 
exists in the periodic Lorentz gas only in the regime of very high densities, where a simple Boltzmann approximation 
may not necessarily be expected to work. We first remark that in a Lorentz gas there are no density corrections due to 
screening of particles, or to many-particle collisions, as they are contained in an Enskog equation p9|. However, the 
existence of an offset at large densities may be related to the fact that the Boltzmann approximation represents only 



the first term in a series expansion in n, which is known as the density expansion of kinetic theory. For the random 
Lorentz gas, such a density expansion has been carried out exphcitly in Refs. jsj,^. It has been found that there 
exist higher-order terms being logarithmic in the density corresponding to so-called ring collisions, which diminish 
the Boltzmann diffusion coefficient quantitatively. However, these logarithmic corrections are difficult to see in the 
functional form of the diffusion coefficient. It is not known to us whether a density expansion has been performed for 
the periodic Lorentz gas. But the offset we find might be associated to the existence of such higher-order corrections 
in the density, and the apparent Boltzmann-like linear behavior may be related to the observation that in similar 
systems higher-order corrections do not change the functional form of the diffusion coefhcient in a drastic way. 

Because of the functional relation between number density n and gap size w Eq. (ffl) , and according to the Boltzmann 
approximation Eq. (||), there should exist a respective approximate linear dependence of the diffusion coefficient D on 
w, as can in fact be seen in Fig. 0. The quadratic term in w which appears in Eq. {M turns out to be quantitatively 
at most 6% for largest w, and, after a proper adjustment, cannot be seen as any deviation to a linearity in w. 

V. FINE STRUCTURE OF THE DIFFUSION COEFFICIENT 

To learn more about the very detailed dependence of D on the density n, or alternatively on the gap size w, we 
performed further computer simulations in the region of large w. The results are shown in Fig. pi Here, the respective 
residuals of D have been plotted, that is, the deviation of the diffusion coefficients from a linear fit in w over the whole 
region shown in the figure. At each value of w ten independent runs of more than 2 x 10^ collisions each have been 
carried out yielding the diffusion coefficients shown by the dots. The squares correspond to the averages over the ten 
runs, where the size of the squares indicates the size of the numerical error. Fig. may be regarded as a signature 
of long-range correlations in microscopic scattering events. A certain subclass of such scattering events are the ring 
collisions mentioned above leading to a series of logarithmic corrections, and divergences, in the density expansion of 
the diffusion coefficient. It might be conjectured that the existence of these divergences in the density expansion of 
kinetic theory is due to the problem of approximating a diffusion coefficient, which is in fact a function as complicated 
as the one shown in Fig. H, in form of a simple series expansion, as has already been remarked in Ref . pO] . 

VI. CONCLUSIONS 

We have numerically computed the density dependence of the diffusion coefficient in the two-dimensional periodic 
Lorentz gas. A comparison of our results with the analytical approximation of Machta and Zwanzig shows that their 
argument is only asymptotically correct in the limit of high densities. We have furthermore discussed how their 
approximation can be corrected systematically by including correlations in the hopping mechanism. 

As an alternative way to understand the structure of the density-dependent diffusion coefficient, we compared our 
numerical results to a simple Boltzmann approximation. We found that on a coarse scale there exists a dynamical 
crossover from a function linear in 1/n for lower densities, as qualitatively be described by the Boltzmann diffusion 
coefficient, to a behavior at high densities which can qualitatively be understood by the Machta-Zwanzig argument. 
Our two different points of view of understanding the density-dependent diffusion coefficient - namely, on the basis of 
the Machta-Zwanzig argument, and by a Boltzmann approximation -, are not at all contradictory. The complicated 
microscopic scattering processes responsible for the failure of the Machta-Zwanzig approximation apparently just 
superpose in a way such that the resulting diffusion coefficient is approximately linear in 1/n on a coarse scale. A 
similar relation between nonlinear corrections and a linear response has been discussed in Refs. p7| , |3q | with respect 
to the existence of Ohm's law in a periodic Lorentz gas with an external field. That this picture is too simple is 
demonstrated by the fact that on a very fine scale we find deviations from such a simple linear behavior. These tiny 
fluctuations may be regarded as the signature of specific microscopic correlated scattering processes in the deterministic 
diffusion coefficient. 

A dynamical crossover between different asymptotic laws for a parameter-dependent diffusion coefficient has already 
been found and discussed in related periodic one-dimensional chaotic maps |10| , |ll| . As in our discussion above, for 
these maps two different random walk models have been used to describe the corresponding asymptotic behavior: For 
small parameters, the respective random walk basically depends on the hopping probability of a particle leaving a 
cell, thus corresponding to the Machta-Zwanzig approximation, whereas for larger parameters the respective random 
walk scales with the distance a particle travels per time step, thus corresponding to the Boltzmann approximation, 
where diffusion is proportional to the collision length of the moving particle. We furthermore note that the diffusion 
coefRcients in this kind of maps have been found to be fractal with respect to variation of the control parameter, and 
that specific correlated microscopic scattering events could be identified as being responsible for this fractal structure 



L0[ . The same phenomena have been encountered in related two-dimensional multi-baker maps, which are simple 
toy models of the periodic Lorentz gas [|l6| . 

Whether the density-dependent diffusion coefficient of the periodic Lorentz gas is in fact a non-differentiable func- 
tion of the parameter, as it is the case for the one-dimensional maps discussed above, remains an open problem. It 
seems to us that this question cannot be answered by performing straightforward computer simulations, but that, 
instead, more refined techniques to compute the diffusion coefficient are necessary, as they have been developed, for 
example, in Ref. PQ | in case of one-dimensional maps. A solution to this problem appears to be important for a 
complete understanding of deterministic diffusion in the periodic Lorentz gas. 
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APPENDIX A: PROBABILITY FOR A COLLISIONLESS FLIGHT ACROSS A TRAP 

In this part of the Appendix we describe a simple analytical approximation by which the probability pd for a 
coUisionless flight across a trap can be computed. The result is shown in Fig. ^ (a). Our approximation is based on 
straightforward applying the Machta-Zwanzig phase space argument described in Section II. It says that 

where A is the area of the trap, 27r is the measure of the velocity space of the trap, S is the region of the boundary 
of the trap across which coUisionless flights are possible, t; = 1 is the velocity of the moving particle, r is the average 
time at which a particle leaves the trap, and 

= djTi-v Q{"f) (A2) 

is the average flux of particles with velocity v_ parallel to the normal n at one entrance, with 7 being the angle between 
n and t;, and ^(7) = cos 7 being the invariant probability distribution of 7 in equilibrium. Here we have already 
assumed that, in an approximation, (j) is completely independent of its position at the boundary. In Ref. pl| it has 
been noted that 

r = ^ , (A3) 

6w 



where w is the size of an entrance, such that Eq. (Al) reads 



Pc. = f^ . (A4) 

Thus, it remains to compute S and (j). The relevant variables for this purpose, and the geometry, are defined in Fig. 
g, where, because of symmetry, we have reduced the problem to one entrance. In a second approximation, we assume 
that regions of size e, as the one indicated in the figure, are the only parts of the boundary which contribute to free 
flights, e is precisely defined as the section on the entrance line between the two horizontal lines which touch the 
upper scatterer tangentially, and which go from the lower edge of the left entrance to the lower edge of the right 
entrance, respectively. Because of symmetry, there are six regions like this for each trap leading to S" = 6e. e can 
simply be calculated to 

e ^ ^(^ +;j 1^-1 . (A5) 

2sin(7r/3) sin(7r/3) ^ ' 

Note that only for e > coUisionless flights across a trap are possible. In other words, e = or, by employing Eq. 

{Ml 



^^^2(sin(V3)+l)_2^^^347 (A6) 

v3 

marks just the onset of the existence of coUisionless flights with respect to varying w. 

It now remains to compute the flux (f> in these regions. To do this, we need to know the hmits of the integral in Eq. 



(A2). As a very crude approximation, wc assume that these boundaries are determined by the angles /3 and S such 

^7 cos 7 



that 



= (5 + 1/2 cos(7r/12) sm{2S) . (A7) 

Here, /3 = tt/6 is identical to the angle between a horizontal velocity vector v and the normal n of the entrance. 6 
gives a bound for the directions the velocity v_ can take such that a particle going in through the region of size e can 
perform a coUisionless flight across the trap. By a straightforward calculation, it can be found that 

2ecos(7r/6) 

d = arctan- — ; — -— ; — ; — 7-— . (A8) 

2 + u;-2(e+l)cos(7r/3) + esin(7r/6) ^ ' 

In summary, we have computed the probability pd for coUisionless flights across a trap as 

Pcf = ^ , (A9) 

w 



where e is given by Eq. (Ka), and is given by Eqs. (A7) and (A8) 



APPENDIX B: HIGHER-ORDER BACKSCATTERING APPROXIMATIONS OF THE DIFFUSION 

COEFFICIENT 

Here we give the explicit formulas for higher-order approximations of the diffusion coefficient in the spirit of Eq. 
(|6|), where this equation represents already the first order correction. The basic idea of these corrections is discussed 
in Section 



IHC 



The underlying picture for all these approximations is the Machta-Zwanzig assumption of diffusion as a hopping 
process with frequency r^^ on a hexagonal lattice of traps, where the lattice sites are separated by a distance of I. 



We use in the following the symbolic dynamics introduced in Section HI C to denote the transition probabilities of a 
particle which hops from one initial trap to some other trap on the lattice within a certain time interval of multiples 
of T. For the particular time interval of 3t, wc find that the particle can move 

- with a probability pi = p{zz) + p{zl) + p{zr) + p{lz) + p{rz) over a distance li = I = {2 + u;) v3 

- with a probability p^ = p{lr) + p{rl) over a distance I3 = \/ll 

- with a probability py — p{ll) + p{rr) over a distance ^3/ = 21. 

The corresponding second-order diffusion coefficient then reads 

D2^il2T)-'[pill+P3ll+P34,] . (Bl) 

In the same way the third-order approximation is obtained. Here, we have a collection of 27 probabilities correspond- 
ing to the different lattice sites where the particle can move within a time interval of 4t. These transition probabilities 
are labeled by all possible combinations of the three symbols 1, r, and z corresponding to sequences of three symbols. 
Thus, the particle can move 

- with a probability po — p{lzz) + p{rzz) + p{zlz) -f p{zrz) -f p{zzz) over a distance ^0 = 

- with a probability p2 = p(lll) +p{llz) +p{lrz) +p{lzl) +p{lzr) +p{rlz) +p{rrr) +p{rrz) +p{rzl) +p{rzr) +p{zU) + 
p{zlr) + p{zrl) + p{zrr) +p{zzl) +p(zzr) over a distance I2 = ySl 

- with a probability pn — p{lrl) -f p{rlr) over a distance /41 — 3/ 

- with a probability P42 — p{llr) + pp{lrr) + p{rll) + p{rrl) over a distance Z42 = (1.5 + \/3)l. 

The corresponding second-order diffusion coefficient then reads 



i?3 = (16r)-i[p2?2 +m'4i +PA2II2] ■ (B2) 

Eq. (Bl) thus corresponds to the second-order approximation plotted in Fig. |g (a), and Eq. (B2) yields the respective 
third-order approximation. 
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FIG. 1. Geometry of the periodic Lorentz gas. A point particle scatters elastically at hard disks arranged on a triangular 
lattice. The cross-hatching indicates a single triangular trapping region. The smallest interdisk distance is denoted by w. 
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FIG. 2. Diffusion coefficient D as a function of the gap size w. The dashed line represents the Machta-Zwanzig random walk 
model Eq. (|2|), the different symbols refer to single data points obtained in the literature from various methods (see text). The 
crosses connected with lines are our new results from computer simulations. All of our data points have error bars smaller than 
the symbols. Error bars of the other data points have been included as far as available. 
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FIG. 3. Correction of the Machta-Zwanzig approximation by coUisionless flights: (a) probabil ity Pcf of coUisionless flights 
across a trap, numerical results (dotted line with crosses), and the analytical approximation Eq. (A9) (bold line) (b) diffusion 
coefficient D, numerical results (bold line with crosses) in comparison to the Machta-Zwanzig approximation Eq. (H) (dotted 
line), and compared to the correction of Eq. (^ by including pcf (dashed line). 
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FIG. 4. Correction of the Machta-Zwanzig approximation by backscattering: (a) backscattering probability pbs as a function 
of the gap size w. For larger w the single data points are plotted by symbols and are connected with lines, for smaller w 
only the lines are shown. The dotted line corresponds to the value of 1/3 of equal probability for any gap, as it is assumed 
in the Machta-Zwanzig approximation. The inset is a half-logarithmic blowup of the initial region for small w. The bars 
included in the figure refer approximately to the regions of different slope in the main figure, (b) numerical results for the 
diffusion coefficient (bold line with crosses) in comparison to the Machta-Zwanzig approximation Eq. (0) (dotted line), and to 
the correction via pbs Eq. (W) (dashed line). 
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2x/w 

FIG. 5. Initial conditions yielding backscattering for w — 0.1. Each dot represents an initial condition for which the particle 
leaves the trap through the same aperture where it entered. 2x/w is the position of a particle on the line of the entrance, and 
a is the angle of flight of an in-going particle to an axis perpendicular to the entrance. 
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FIG. 6. Systematic approximations of the exact diffusion coefficient in higher order: (a) combined approach of employing the 
numerical probabilities of coUisionless flights and of backscattering within a hierarchy of analytical expressions of the diffusion 
coefficient. The flrst-order approximation is given by Eq. (n) by using the probabilities pd and pha of Figs. H (a) and H (a). The 
second order corresponds to Eq. (Bl), and the third order corresponds to Eq. (B2|). They include higher-order backscattering 
events, as described in the text and in Appendix H. (b) diffusion coefficients obtained from lattice gas simulation by using 
probabilities of correlated hopping sequences as calculated in Lorentz gas simulations (see text). 
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FIG. 7. Diffusion coefficient in the periodic Lorentz gas with respect to the inverse of the number density of the scatterers 
n. The Boltzmann approximation Eq. (H) is shown after subtracting an offset of AD — 0.025 (thin bold line). It is compared 
to the numerically exact results (thick bold line with crosses), and to the Machta-Zwanzig approximation Eq. (0). 
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FIG. 8. Residuals of D), that is, the difference between the numerically computed diffusion coefficient D and a linear fit over 
the range of the data presented in the figure. The dots indicate the results of ten independent runs for each value of w, the 
squares represent the algebraic averages over these ten runs. The size of the squares is equal to the magnitude of the numerical 
error. 
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FIG. 9. Geometry to compute the probability of coUisionless flights Pd across a trap, e labels the size of one of the two 
regions at any of the three entrances of the trap where ingoing particles can perform coUisionless flights. 5 is related to an 
approximate bound for the directions the velocity of such a particle can take such that it performs a coUisionless flight. 
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